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The Navier-Stokes equation describes the deterministic evolution of incompressible fluids. 
The effects of random initial conditions on solutions of this equation are studied. It is 
shown that there is an infrared stable fixed point accessible within the epsilon expansion, 
about d c = 4, with a particular choice of viscosity. The behaviour of the usual viscous 
fluid is obtained as a relevant perturbation about this fixed point, with the usual viscosity 
serving as T — T c . 
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The equations of fluid mechanics are deterministic. The motion of an incompressible 
fluid, in the absence of viscous forces, is governed by 

1 

d t u + u-Vu = — Vp, (1) 
P 

and the incompressibility condition, V • u = 0. Here u is the velocity of a fluid element, 
p is the pressure, and p is the constant density of the fluid. In fact, the pressure term is 
determined up to a harmonic function by the velocity field using V-w = 0. This equation has 
no obvious expansion parameter, since all the terms, with precisely determined coefficients, 
are determined by Newton's equation, F = ma — as a corollary, the equation possesses 
Galilean invariance. 

In principle, all details of the fluid's motion are determined by completely defined 
initial conditions. The ideal aim of a theory of turbulent motion would be to establish a 
statistical description analogous to the kinetic theory of gases, based on the deterministic 
dynamics of fluids. The statistical description may be thought of as an ensemble average[l], 
over different experiments corresponding to different initial conditions. It is the aim of this 
paper to perform such an average over random initial conditions for incompressible fluids. 

The approach taken here is analogous to the standard infrared stable non-trivial fixed 
point found in the e expansion for </> 4 theory, about d c = 4 [5]. In that case, one considers the 
massless theory about d c = 4, expands in e = d c — d, and one finds a fixed point with e oc g. 
Correlation functions at T > T c , in the critical domain, i.e. with the correlation length 
large compared to the cutoff, but still finite, exhibit strong scaling. Physics at length 
scales small compared to the correlation length, but large compared to the underlying 
cutoff, can be expressed in terms of universal correlation functions in which the unit of 
length is provided by the correlation length. 

Turbulence, in general, involves the transfer of energy from large-scale mean motions 
to viscous dissipation at small scales. The case of isotropic homogeneous freely decaying 
turbulence considered in this paper, however, involves no energy pro duct ion [1], and hence 
one need not take into account the so-called 'energy cascade', only viscous dissipation at 
small scales. Ultimately, all turbulence dies away due to non-universal dissipative viscous 
forces, so the theory of interest is only expected to apply in the inertial regime where 
the viscous forces dominant at long distances can be neglected. These dominant viscous 
forces determine the length scale of decay, so in the inertial range we are interested in 
fluid velocity correlations at distances long compared to the underlying molecular cutoff, 
but short compared to this viscous length scale. Compare this description of the inertial 
range with the description given above of (j) 4 theory with T — T c small but non-zero. In 
the inertial range then, the viscous forces of interest are higher terms in the expansion of 
viscous forces in powers of wavenumbers (and, in principle, nonlinear terms in the velocity 
itself — these will turn out to be irrelevant) . We will find that there is a critical point in this 
'hyperviscous' theory, with the viscous theory in the inertial range described by universal 
correlation functions in which the unit of length is provided by the standard viscosity. In 
this sense, our approach is modelled closely on the successful e expansion framework of 
critical phenomena. 

Now, while the operational meaning of the ensemble average of physical interest was 
described above, previous theoretical investigations in a field-theoretic framework have 
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added a 'random force' term to the basic equation — the resulting exponents of turbulent 
flow depend strongly on the scaling form assumed for the force-force correlation function, 
contradicting the observed universality of the exponents [4]. By the fluctuation dissipation 
theorem, the amplitude of the random force has to be related to the viscous forces for 
thermal equilibrium to be reached at long time scales. This correlation implies that the di- 
mensionality of the order parameter (in this case the velocity field) depends on the random 
force. This is an explanation for the non-universality of theoretical computations [2]. 

Several years ago, Liao[2] examined the effects of averaging over all solutions of eq. (1) 
in a field theoretic setting, and obtained significant results on the spectra of fully developed 
isotropic turbulence. My analysis is similar in spirit to Liao's approach, but is based on 
Cardy's calculational framework for the effects of random initial conditions on deterministic 
diffusion equations with nonlinear terms [3]. 

The initial conditions I consider are the simplest allowed by locality and incompress- 
ibility: 

(u l (x, 0)u j (y, 0)) = -g(AS ij - d l d J )5{x - y), 

with A = V 2 . It is important to keep in mind that there is an underlying cutoff in the 
problem, even though our computations will be done with dimensional regularization. I 
shall assume (u l (x, 0)) = 0. From these initial conditions, and eq. (1), observe the following 
engineering dimensions: [u] = Ek~ x , [g] = [u\ 2 k~^ d+2S) . In an isotropic incompressible 
homogeneous fluid, the viscous forces are the only non-universal terms in the equation. I 
am interested in the inertial regime, where the viscous forces are negligible, yet the form 
of the viscous terms determines the dimension of [u]. Assume that the viscous force is of 
the form upA^u + lower order, Then [up] = Ek~ 2 ^ . In principle, one could eschew the 
analyticity in momenta assumed implicitly here, but the framework of renormalizability 
must then be considered ab initio. Following Cardy[3], I expect that the only quantity 
that requires renormalization will be g, since the evolution equation is deterministic, and 
the only functional average performed is over the initial conditions. I will henceforth set 
i/0 = l. Then 

[g] = fc4(0-l)-d 

so for d c = 4, (3 = 2. As explained above, this is what we expect in order to find the critical 
point suited to a description of the inertial regime: The standard viscous force, vAu, will 
appear in our theory as a 'massive' deformation away from the critical theory, with the 
hyperviscous force, z/ 2 A 2, u, determining the critical theory. Of course, a standard mass 
term would violate Galilean invariance for a fluid, so in fact this analogy is quite precise. 

As mentioned above, the compatibility of eq. (1) with the incompressibility conditions 
requires that p is related to u by 

--Ap = diUjdjUi, 
P 

in order that initial conditions satisfying V • u(x,t = 0) = produce a flow satisfying 
V • u(x,t) = 0. This expression is dimensionally correctfl]. (Compare with Landau and 
Lifshitz[6, eq. (31.4)].) Only the gradient of the pressure enters eq. (1), so there is no 
ambiguity in inverting the Laplacian for our purposes. The initial conditions are explicitly 
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local, the inversion of the Laplacian needed for the pressure term is merely a consequence 
of the infinite velocity of sound in an incompressible fluid. 

Treat the nonlinear terms in eq. (1) formally as perturbations. This leads to an 
iterative solution to the equation which has the structure of a tree graph, as is standard [5]. 
The trunk of the tree represents the solution at time t, u(x,t), and the time-ordered 
branches all terminate at u(x{, 0). The ensemble average of interest has the effect of closing 
the branches into loops, with attendant factors of g. It is then clear that g is the expansion 
parameter of interest — remarkably, the ensemble average has a perturbative expansion, 
even though the actual equation has no such perturbative parameter. 

Define U d (k) = 5 lJ — k % W jk 2 . The formal solution of the equation then has the 
following structure: 



u\k,t) =e- tk4 u l (k,0)-iU il (k)k J [ dt'dp 

Jo 



1(! (t-t')k* e -t'ip*+(k-pr)„i 



u L (p : 0)u J (k-p, 0) + . 



(2) 

Any correlation function of the form (f| u ki (xi, ti)) can be evaluated in a straightforward 
fashion, using the assumed form of the initial conditions. To find a fixed point, the renor- 
malization of g must be computed. It is elementary to observe that at order g p the degree 
of divergence of an n point function is 



5 = d - (d - 3)n + (d - 4)p. 

It would appear then that there are three possible primitively divergent functions at d = 4 
(since the one point function is trivial by Galilean invariance). However, this is not actually 
the case, due to the incompressibility condition. This condition has a drastic effect on the 
vertices, resulting in explicit factors of momentum coming out of loop integrals, apparent in 
eq. (2). The actual degree of divergence of the two point function is only logarithmic, and 
the other functions are actually finite. The divergent function can be rendered finite with 
just a renormalization of g, due to this fact. A corollary is that the case of a compressible 
fluid is not a trivial extension of the present approach. 

The integrals that arise in loop calculations for arbitrary times are unfortunately not 
easily performed, due to the quartic character of the hyperviscous force. To compute 
renormalization group functions, I need only the explicit value at E = 0, going to an 
energy-momentum representation. The divergent contribution to the two point function 
at 0(g 2 ) comes from the integral 



1 = 



d d p p 2 (q-p) 2 U Lr (p)U ms (q-p) 



(2tt) c 



(p* + {q _ p) 4y 



There are other contributions at 0(g 2 ) that are manifestly finite, and difficult to evaluate. 
The two point function (excluding uninteresting finite parts at 0(g 2 )) is (e = 4 — d) 



(u i (q,E = 0)u\-q,E = 0)) 



2 ng(g) / fi- 
= 167T g t — 1 + Cg — 

q 2 =M 2 q \ e 
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I have rescaled g to absorb factors of 167r 2 , and C = 7/24. Define the dimensionless 
renormalized coupling by g r = g^~ e (l + Cg^-^j , then 

d 

^~d\x 9r ^ ^ 9r ^ = ~ e9r ~ Cg *~ 

There is a zero of the beta function at g* = — e/C, with /3'(g*) = e > 0, implying infrared 
stability. The sign of g* is important: The fact that it is negative signifies that the critical 
point is associated with a multi-modal distribution, which is physically quite appropriate 
for turbulence — we do not observe turbulent behaviour starting spontaneously in beakers 
of water from thermal fluctuations, an example of a unimodal set of initial conditions. 

Various simple scaling properties at the critical point are easily derived in the standard 
renormalization group framework[3]. For example, at fixed q, E, 

(u(q,E)u(-q,E))^^ g ;, 

with 7 = {(3 + eg r )/g r , and 7* = 7(0*) = e. It follows that 

(u(x, t)u(x, t)) ~ t~ , 

independent of e. Similarly, one deduces that 

(u( Xl ,t) . ..u(x n ,t)) = t- 3n / 4 G^(\x % - ^-It" 1 / 4 ), 

which is essentially just dimensional analysis taking into account the fact that u is not 
renormalized. Notice that in principle we could also have found universal exponents by 
considering the response functions corresponding to variations of the initial conditions, as 
studied in Ref. 3. In the present case, the incompressibility condition implies that these 
functions are finite, hence do not exhibit any nontrivial scaling. The study of these response 
functions may lead to a better understanding of the convergence of the e expansion in this 
theory. 

The physical case of interest is really the behaviour of the theory away from the critical 
point, perturbed by the usual viscosity term. As is quite standard[5] the behaviour of 
correlation functions at distance scales small compared to the finite correlation length, but 
long compared to the cutoff, is obtained from the critical theory with appropriate resummed 
operator insertions. This analysis poses no new conceptual problems in the present theory, 
but is computationally rather involved due to the quartic energy denominators. However, 
even if the integrals are analytically intractable, the numerical solution at e = 1 should be 
quite feasible. 

In summary, I have presented here concrete evidence within the standard e expansion 
for a nontrivial infrared stable fixed point governing, in the inertial range of interest, 
statistical properties of solutions of incompressible fluids when the initial conditions are 
multi-modal, without the introduction of any random force. Since numerical techniques 
to cope with the e expansion are very well studied, it should be possible to derive detailed 
predictions for experimental observations of universal behaviour in turbulent flows starting 
from correlation functions at this fixed point. These issues, and especially the extension 
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to compressible fluids, will be discussed elsewhere. Numerical simulations of fluids could 
be directly used to test the hyperviscous fixed point at criticality. 
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